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Random Notes on Matrices 

K. Goldberg 

Three matrix problems are considered in this paper: Bessel functions as limits of de- 
terminants, finding all optimal strategies of a matrix game with nonzero value, and con- 
ditions for matrices to have equal principal minors. 

The three sections of this paper have little relation except that they each concern matrices 
and each was suggested by the work of a colleague at the National Bureau of Standards. 

In the first section we prove that the limit of a certain sequence of determinants is the 
modified Bessel function A s (x). 

In the second section we show that a method for computing a pair of optimal strategies 
for certain square matrix games yields the unique optimal strategy for one player if the optimal 
strategy for the other player is totally positive. This can result in an appreciable saving in 
the work of computing all the optimal strategies of an arbitrary matrix game with nonzero 
value. 

The third section contains an amplification of a result which appeared in this journal as 
well as a necessary and sufficient condition that two matrices have equal corresponding princi- 
pal minors. 

1. Bessel functions as limit* of determinants. Let (a ) be an infinite, triple diagonal, symmetric 
matrix with 

< i -i=a 11 -an=l-(2 5 +l)i- 1 +i- 2 ^- 1 ) (1.1) 

for some polynomial p(t) and some real number s which is not a negative integer. Let D n (t) be 
the characteristic polynomial of the first n rows and columns of (dy). We shall prove 

Theorem 1. lim D n (a n + l— n~ 2 z 2 ) = A s (z), I^K 00 , 

where A s (z) is the Jahnke-Emde modified Bessel function with expansion 

oo / -i\kr?2k 

As( ' 3)=: § A k k\(s+k)(s+k-l) ... (s+1)' N<°°- 

The theorem was suggested in considering a conjecture of J. Todd. 

Since (a 7j ) is triple diagonal and symmetric we can use the first equation in (1.1) to obtain 
a simple recursion for D n (t). We have D (t) = l, D x (t)=t—a n and 

D n (t) = (t-a nn )D n ^(t)-(a n -a nn )D n . 2 (t) n=2, 3, . . . . (1.2) 

Thus 

Z> n (a u + l)=l. 

We shall expand D n (a n J r l~-n~ 2 z 2 ) around a n + l, and to this end we define 

Z%"=Wan+l). 

By repeated differentiation of equation (1.2) we get the following recursion for 
D^:D^=D^ ) = \ ) D^ k) =0 if k>n and 

D<*>-2D<*2 ,+Di^+n-' (2s+ 1) (P™i-D™J -n^n' 1 ) (D^-D^) =kD^ 

for &= 1,2 . . . ,n and ft — 2,3, .... This recursion can be solved for D ( n k) yielding 

D\ k) = n 2k A~ k II (j+s) - 1 +n 2k - l J k +n 2k - 2 q k (n- 1 ) 
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for all n and k<n, some polynomial g k (t) of degree 2k— 2, and 

a =d/2Jfe\ r(fc+H-i)r(H-i) * , _ 2 / 2 A r(2j+2s-i) (2j+i)(2H-i)+2 P (0) 

a * *U/ r (2AH-2*+l) U \jjTU+s)T(j+s+l) 2j-l 

for all k. All that we need to know about d k is that it is bounded and this is clear from the most 
crude inequalities. Let d be the bound: 

\d k \<d, 4=1,2, .... 
Denote the partial sums of A s (z) by 

(— l)*g* 



a s „0)=]C- 



"fej4*-i!(s+fc)(«+fc-l) . . . (*+D 
Then we have 

^(an+l-^- 2 3 2 )-A sw (s)=g^ ^■#~4*( 5+ jfc)( g+ jk-l) . . . (s+l)/ 

so that 

Z?»(aii + l-w-V)-A, n (g)<n- 1 de , * ,a + o(l), |z|<>. 

Since A sw (^) tends to A s (^) for all |^|<C °° as w gets large the theorem is proved. 
Note that 

A_^(s) = cos z } 

Ay 2 (z) = z~ 1 sin 2, 
and in general 

A s (2) = r( s +i)(^/2)- s J s (2), 

where J* (2) is the Bessel function of order s. Thus our theorem can be interpreted as meaning 
that if z } ,z 2 , . . . ,z n are the n smallest roots of J s (z) then the roots of D n (f) are approximated 

by 

dn+1—n 2 z\ i=l,2, . . . ,n, 

with the best approximations being for the smaller z it 

2. On finding all optimal strategies of a matrix game with nonzero value. Let B be a payoff 
matrix with nonzero value. Let K x and K 2 be the extremal sets of the optimal strategy spaces 
of the row player and the column player, respectively. 
A result due to L. S. Shapley and K. N. Snow [1] states 

Theorem. Let x and y be optimal strategies for the row player and the column player 
respectively. A necessary and sufficient condition that xeK x and yeK 2 is that there exist a 
nonsingular submatrix A of B of order n such that, if e is the 1 by n vector of all ones, 

v=[eA~ l e T ]- 1 

x=veA~ l (2.1) 

y=ve(A T )~ 1 

where v is the value of the game, x is the vector obtained from x by deleting the elements 
corresponding to the rows deleted to obtain A from B, and y is the vector obtained from y by 
deleting the elements corresponding to the columns deleted to obtain A from B. 

We shall prove that if either x\> (that is every element in x is positive) or y^> then the 
other represents the unique optimal strategy for A for the corresponding player. In other 
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words if, say, x^>0 then no element of K 2 other than y can be obtained from any submatrix 
of A 

This result can be used to reduce the number of submatrices of B that have to be inspected 
to exhaust K t or K 2 . 

Our result takes the following form : 

Theorem 2. Let A be a nonsingular payoif matrix of order n such that the nonzero row and 
column sums of A~ l are all of the same sign. If none of the row (column) sums of A~ l are 
zero then the row (column) player has a unique optimal strategy. 

We shall say that 1 by n vectors x and y represent optimal strategies for the row player 
and the column player of A respectively if 

£o>0, #o>0, x e T =y Q e T =l 
and for some number v , the value of the game, 

XoA>v e>y A T . 
Let v, x and y be as in (2.1). Then xe T =ye T =l and 

xA=ve=yA T . 

By hypothesis x>0 and y>0 so that x and y represent optimal strategies and v is the value 
of the game. 

We shall assume that x\> and prove that y represents the unique optimal strategy for the 
column player. The proof for the other case follows in exactly the same way so this will prove 
Theorem 2. 

Suppose y is a non-negative 1 by n vector such that y Q e T — 1 and 

ve>y A T . (2.2) 

Since £>0 we can find a 1 by n vector u such that the matrix 

X=x T u-(A T )~ 1 
is non-negative. Then (2.2) implies 

veX>y A T X. 
Computing these two vectors we get 

veX= vex T u— ve (A T ) ~ l = vu— y 
and 

y A T X=y A T x T u—y Q A T (A T )- 1 ==vyoe T u—y =vu—y Q . 

Therefore we have y$>y. Since the sum of the elements of these two vectors are equal we 
must have yo=y, as desired. 

3. Matrices with equal principal minors. In a previous paper in this journal [2] we stated a 
result which is a special case of the following 

Theorem 3. Let {a tj ) be a matrix of order n with complex elements and the properties 

dijdji is real and non-negative, i,j=l,2, . . . ,n, (3.1) 

_ , (3.2) 

a hu a hH ' ' • a hh — (li 2 i i a hi2 ' ' ' a h*k 

for all fc = l,2, . . . ,n and i j =l,2 J . . . ,n. Then (a w ) has the same principle minors as the 
Hermitian matrix (b tj ) with 

b tj = (sgn^dij) {difiji) X (3,3) 

and so has only real characteristic roots. 
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The symbol x means the complex conjugate of x, and the symbol x 1/2 means that square root 
of x with non-negative real part. 

We shall prove this theorem by showing that 

for all Ar=l,2 7 . . . ,n and i j =l,2 J . . . ,n and the result will follow as in the previous paper. 
From (3.1) we see that either a t ja ji =0 or 

Uji — Cij&ij 

for some positive real c 4i . These c^ satisfy 

Cij—Cjt 

and, from equation (3.2), 

c hh c hH - - - Ci kh = ^ (3. 5) 

whenever a ili2 a i2iz . . . a ijcil does not vanish. 
Therefore, either b ij =b ji =0 or, 

so that 

as desired to prove that (6 y ) is Hermitian. 

Choose an arbitrary ordered subset i h i 2 , . . . ,i k of {1,2, . . . ,n}. Let 6 = 6^6^3 . . . b ikil 
a=a ili2 a i2is . . . a ijciv a' = a i2h a i?)i2 . . . a ilijc . 

If b vanishes then, by (3.3), either a or a' vanishes. By (3.2) this means that a vanishes, 
so b — a. 

If b does not vanish we have 

by combining equations (3.5) and (3.6). 

Therefore b = a in general and we have proved equation (3.3) and so the theorem. 

We are naturally led from this theorem to ask for a necessary and sufficient condition that 
two matrices have the same principal minors. We shall find such a condition based on a little 
known expansion of the determinant of a matrix. 

Let A=(a tj ) be a matrix of order n. Choose any set of k distinct indices [i u i 2 , . . . ,4} 
with l<ij<n. Then define the symbols 

itk=l, and (3.7) 

\HMj • • • )^k\ A == 2-i a iii ff( i ) a i (r{ iy (T 2 (1) • • • a i <r k-i {1) i 1 
a 

if fc=2,3, . . . ,11, where the sum is taken over all full cycles a on {1,2, . . . ,fc}. 
Given an arbitrary partition of {1,2, . . . ,n} 

T== {^11^12? • • • ^l*i ){ ^21^22? • • • Mk 2 J - • • {W>to2> • • • j'W^} 

define the symbol 

TA={—l) m {huh2, . . . M* x ) a^iM) • • • Mk 2 } A • . . {imuim2, . . . Mk m } A . (3.8) 
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Then the determinant of A may be expanded in 

det^=(-l)*2>,i (3.9) 

where the sum is taken over all partitions r of {1, 2, . . . , n}. 

To prove that this expansion is valid we need only show thai the generic term in the usual 
expansion of the determinant appears in this expansion, that whenever it does appear it docs 
so with the proper sign, and that there are exactly n\ such terms. 

The generic term in the expansion of the determinant can be denoted by 

CL=( l) P &lp(l)&2p(2) • • • Q>np(n) 

where p is a permutation of degree n and (— l) p is +1 or — 1 as p is even or odd respectively. 
Write p in its unique decomposition into cycles 

P— (?11?12 • • • ^1*^(^21^22 • • • r hk 2 ) • • • (W^m2 • • • ^rnk m )' 

Then the only r A in which a can appear is the one described in (3.8). Since (-l)p=(-l) w-w the 
sign is always right. Finally it is clear that a appears at least once in r,. 

Now, consider the partition ki J \-k 2 -\-. • .-\-k m of n. Let s t be the Dumber of kj equal to i. 
Then there are 

n! 

s x \s 2 \ . . . s n \ki\k 2 \ . . . k m \ 

partitions r of {1,2, . . . ,n} with subsets of order k h k 2 , . . . , k m . By (3.7) a subset with 
ki elements contributes (k f —\)l summands to r 4 . Thus the total number of summands on the 
r.h.s. of equation (3.9) is 



Zj si\s 2 l . . . s n \ h\k 2 \ . . . k m \ Zj sM • • • s n \l 8 



l 2 S2 . . . n Sn 

as desired, the sums being taken over all partitions of n. 
Suppose B=(bij) is a matrix of order n and 

{i h i 2 , . . . ,i k } A ={h,i2, • . • , h}b (3.10) 

for all subsets {i v i 2 , . . ., i k ) of {1,2, . . ., n] . By (3.8) and the expansion (3.9) this implies 
that the corresponding principle minors of A and B are equal. We shall prove the converse. 
Suppose A and B have equal corresponding principal minors. Then (3.10) holds for k=l. 
Suppose it holds for fc=l,2, . . ., m-1. Consider the expansion as in (3.9) of the minor with 
row and column indices i u i 2 , . . ., i n . For every t except r / ={i 1 ,i 2 , . . ., i m } we are summing 
products of equal sums {i[, % 2} . . . , i'k} a= {ii,i 2 , ■ • •, vIc}b with k<m-l. Therefore the equality 
of the minors implies ta=tb which completes the proof by induction. We have proved 

Theorem 4. Let A and B be mat rices of order n. They have equal corresponding principal 
minors if and only if equation (3.10) holds for all subsets [ii,i 2 , . . ., ik} of {1,2, . . ., n}. 

We remark that theorem 4 (and thus theorem 3 as well) also follows from a result uoted 
by A. Ostrowski (footnote 13 in [3]), which is closely related to the expansion (3.8). 
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